MAGNETO-TRANSPORT AND MAGNETO-OPTICAL PROPERTIES OF 
FERROMAGNETIC (III,Mn)V SEMICONDUCTORS: A REVIEW 



Jairo Sinova,^ T. Jungwirth,'^' ^ and John Cerne* 

^ Texas A&M University, Physics Department, 4242 TAMU, College Station, TX 77943-4242 
''University of Texas at Austin, Physics Department, 
1 University Station C1600, Austin TX 78712-0264 
* Institute of Physics ASCR, Cukrovarnickd 10, 162 53 Praha 6, Czech Republic 
^Depatment of Physics, University at Buffalo, SUNY, Buffalo, New York 14260 

(Dated: February 2, 2008) 

Rapid developments in material research of metallic ferromagnetic (III,Mn)V semiconductors 
over the past few years have brought a much better understanding of these complex materials. 
We review here some of the main developments and current understanding of the bulk properties 
of these systems within the metallic regime, focusing principally on the magneto-transport and 
magneto-optical properties. Although several theoretical approaches are reviewed, the bulk of the 
review uses the effective Hamiltonian approach, which has proven useful in describing many of these 
properties namely in (Ga,Mn)As and (In,Mn)As. The model assumes a ferromagnetic coupling 
between Mn d-shell local moments mediated by holes in the semiconductor valence band. 
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I. INTRODUCTION 

The artificial marriage of ferromagnetic and semicon- 
ducting properties is one of the main motivations behind 
the creation and study of diluted magnetic semicondcu- 
tors (DMS) focused on (III,Mn)V materials. This effort 
is being fueled, in large part, by the endless demand for 
technological advances in our society whose global inter- 
dependence in communications and technology is evident 
throughout our daily activities. Moore's law, predicting 
an exponential growth of computational ability, is now 
believed to be in perryl and there is a need for a new 
technological revolution if the fast pace of technology is 
to be mantained. 

Such a technological revolution will take place if three 
major goals can be achieved in the development of fer- 
romagnetic (III,Mn)V materials: i) room temperature 
ferromagnetism, ii) dependence of transport and optical 
properties on the magnetic state, and iii) that the ma- 
terials mantain their fundamental semiconducting char- 
acteristics, i.e. strong sensitivity to doping, external 
fields, and light, and are compatible and integrable with 
conventional (non-magnetic) micro- and opto-electronic 
semiconductor technologies. Since the first observation 
of ferromagnetism at low temperatures in the (Ill,Mn)V 
semicondcutors in the early nineties,^ steady progress has 
been taking place in each of these major fronts. More 
importantly, over the last five years, jumpstarted by the 
discovery of ferromagnetism in (Ga,Mn)As with ferro- 
magnetic temperatures in excess of 100 K,^ the efforts 
in this field have increased dramatically in several coor- 
dinated fronts, both experimental and theoretical. This 
rapid research has yielded many advances, prompting the 
need for an overall review of at least several fast evolving 
parts of this scientific effort. In this article we discuss 
transport and magento-optical properties of these ma- 
terials focusing on the metallic regime and on theoreti- 
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cal approaches based on the effective Hamiltonian model 
which have been succesful in obtaining qualitative and 
quantitative understanding of various transport charac- 
teristics of these systems. Other important areas, such 
as the insulating low-doping regime have been covered by 
other reviews. '^'^ Furthermore, we will narrow our focus 
on the most widely studied and well understood materi- 
als, (Ga,Mn)As and (In,Mn)As, only mentioning briefly 
some of the current prevailing ideas in regards to the 
much lesser understood ones such as the Mn-doped ni- 
trides or phosphides. 

As in any review, a preference must be given to 
the more relevant publications containing the most re- 
cent developments and the ones which have had the 
largest impact in understanding these complex mate- 
rials. Unless the reference section is equal in length 
to the actual review, it is impossible to cite every 
contribution to this field. To remedy omissions that 
may have partly originated from the authors' own work 
perspective, we refer to a more extended and precise 
database of published and pre-print works maintained 
within the "Ferromagnetic Semiconductors Web Project" 
at http://imixl2.fzu.cz/ms. This web page contains a 
large body of theoretical and experimental results (that 
can be obtained in several formats), as well as the most 
extended bibliography database linked to diluted mag- 
netic (ni,V) semicondcutors that we are aware of. We 
encourage the reader in need of a more detailed bibliog- 
raphy to use this resource. 

We also note that a few review articles on other as- 
pects of the physics of DMS, which we will not discuss 
extensivelly in this review, have appeared over the past 
few years and may help the reader who wants to extend 
the narrow scope of this article. Given however the rapid 
developments over the past two years alone, the reader 
should be warned that these reviews, although cxcclcnt, 
may not sufice to give the most current experimental 
developments and theoretical understanding and justi- 
fications. Several extended papers on the physics of fer- 
romagnetism and proposed mechanisms, some of which 
we briefly review here, are shown in Refs. 2,5-7. A 
theoretical review of the widely used phenomenological 
Zcner mean field model, highly succesful in (Ga,Mn)As 
and (In,Mn)As studies, is presented in Ref. 8 and 9. 
A large body of work up to 2001 in density functional 
first principles calculations in (III,V) DMS is also avail- 
able in Ref. 10. Another theoretical approach based 
on lattice models and dynamical mean field treatment of 
disorder, most succesful in the low-doping regime, is re- 
viewed in Ref. 3. The ever important effects of disorder, 
omnipresent in these materials, is reviewed in Ref. 4. 
A general discussion behind the technological motivation 
for DMS has been presented in 11, 12, 13, and 14, 
where optical isulators and magnetoresistance elements 
are explored. 

We organize the rest of this article as follows: In Sec. 

I A we present a short general overview of DMS from the 
experimental facts point of view. In Sec. II we briefly 



discuss the theoretical model used and its experimental 
basis to describe the metallic regime. Sec. Ill reviews 
the transport properties of DMS such as diagonal re- 
sisticity, magneto-resistance anisotropy, and anomalous 
Hall effect. Sec. IV covers the different optical and 
magneto-optical properties of the systems such as optical 
absorption, magnetic-circular diachroism, Raman scat- 
tering, and cyclotron resonance. A short summary and 
outlook with different open questions is given in section 
V. 



A. The basic picture of DMS 

The basic idea behind creating these novel ferromag- 
netic materials is simple: based on the lower valence of 
Mn, substituing Mn in a (III,V) semiconductor for the 
cations can dope the system with holes; beyond a con- 
centration of 1% there are enough induced holes to medi- 
ate a ferromagnetic coupling between the local 5 = 5/2 
magnetic moments provided by the strongly localized 
3d^ electrons of Mn and a ferromagnetic ordered state 
can ensue. This rudamentary but generally correct pic- 
ture of ferromagnetism in many (III,Mn)V materials is 
of course aquired by peacing together the different ex- 
perimental results obtained through different character- 
ization techniques of the materials. The simplicity of 
the model hides within it a cornucopia of physical effects 
present in these materials such as metal-insulator transi- 
tions, carrier mediated ferromagnetism, disorder physics, 
magneto-resistance effects, magneto-optical effects, intri- 
cate coupled magnetization dynamics, post-growth de- 
pendent properties, etc. 

In any doped semiconductor, an understanding of its 
properties must be preceeded by a study of the doping 
impurities. Since, in this case, the Mn impurity dopants 
provide the local magnetic moments responsible for fer- 
romagetism and, at the same time, provide the hole car- 
riers which mediate the exchange coupling between the 
impurities themselves, understanding the nature of the 
impurities and the different states that they induce is 
even more fundamental in order to establish a relevant 
theoretical model. 



B. Impurities in DMS 

Under equilibrioum growth conditions, the solubility 
of Mn in III-V semicondcutor crystals has an upper limit 
of 0.1 %. Beyond this concentration phase separation 
and surface segregation takes place. To circumvent this 
problem, low-temperature molecular beam epitaxy (LT- 
MBE) technique was applied leading to the first succesful 
growth of InAs and GaAs based DMS.^ The most com- 
mon and stable position of Mn in the host semiconductor 
lattice is on the Ga site.^^"^^ 

An important initial question is what is the neutral 
state of the Ga-substitutional Mn impurity. Experimen- 
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tally, through electron paramagnetic resonance (EPR) 
and optical measurements, ^^'^"'^^ for even very few per- 
centage of Mn (lower by an order of magnitude than 
the concentration marking the matal-insulator transi- 
tion) the only impurity level observed is the d5-Mn^+ 
ionized state, indicating that indeed the strong local- 
ized 3d^ electrons with total spin S = 5/2 are a good 
starting point for theoretical models. The correspond- 
ing 3d''' plus a weakly bound hole neutral state is ob- 
served experimentally only within a narrow range of Mn 
concentrations due to a complete compensation in the 
low Mn doping regime and due to the Mott insulator- 
to-metal transition at high doping. Other candidates for 
the Mn neutral impurity such as 3d^ states, important 
for double-exchange models to be applicable, are not ob- 
served experimentally.^^ Furthermore, the itinerant holes 
originating from the Mn acceptor level have also been ob- 
served through photoemission experiments^^'^^ to be as- 
sociated with the ones of the host semiconductor valence 
band, i.e. they have As 4p-charactcr, which provides 
further evidence for the shallow acceptor nature of the 
substitutional Mn impurity in GaAs and InAs. 

The other impurities present because of the non- 
equilibrium growth process, As-antisites and Mn- 
interstitials, compensate the induced holes and there- 
fore reduce the free carrier concentration with respect 
to the substitutional-Mn density. Annc;aling proce- 
dures, at temeratures slightly lower than growth tem- 
perature, have shown a reduction of the ammount of 
compensation. The initial procedure^ has now been 
modified by different groups^**^^®'^^ and the carrier con- 
centration can actually be tuned precisedly through 
resistance-monitored annealing.^^ Two possible mecha- 
nism for the reduction of compensation have been con- 
sidered: either the number of As-antisites is reduced 
through migration of As atoms to their correct sites, 
or the Mn-interstitials migrate to the surface or to the 
substitutional sites. Early theoretical models assumed 
the former scenario, however, a recent key Rutherford 
backscattering experiments^ has shown that it is the re- 
duction of Mn-interstitials and their out-diffusion to the 
surfacc^*^ which is responsible for lowering the ammount 
of compensation. 



C. Other experimental observations in DMS 

Besides the experimental results revieling what is hap- 
pening in the complicated annealing process and the na- 
ture of the impurities, there are other experimental ob- 
servations that provide key clues for undcrstanting the 
origins and fundamental properties of (Ga,Mn)As and 
(In,Mn)As DMS: 

a) The ferromagnetic behavior of DMS materials is 
only observed above a critical doping level of about 1%.^ 
At lower doping the compensation is almost complete and 
the mediating carriers needed to exchange-couple the lo- 
calized Mn moments are not present and therefore no 



ferromagnetic ordering occurs. 

b) The carrier-induced ferromagnetic nature of the 
ordered state has been demonstrated by field-effect 
(In,Mn)As experiments^" where the carrier concentration 
was tuned by a gate and the corresponding critical tem- 
perature was modified accordingly. 

c) Strain effects, due to the lattice mismatch between 
the DMS layer and the substrate, influence the magnetic 
anisotropy of the ferromagnetic state. It has been shown 
that the ferromagnetic easy axis can be along the growth 
direction or in the plane depending on whether the strain 
is tensile or compressive."'''^ Such phenomena arc inter- 
preted in terms of well known strain effects in the spin- 
orbit coupled valence bands of the host semiconductor. 

c) Both magneto-optical effects, e.g. magnetic cir- 
cular diachroism (MCD),"^^ and resistivity measure- 
ments above the critical temperature, '^^ indicate an anti- 
ferromagnetic coupling between the local 3d^ electrons 
and the valence hole. This so-called kinetic-exchange 
interactions originates from the hybridization between 
the Mn 3d-orbitals and the neighboring As 4p-orbitals 
and is much larger than the pure Coulomb exchange 
interaction.'^^ 

e) In transport measurements, a large anomalous-Hall 
effect (AHE) (see Sec. HID) completely dominates the 
low- field off-diagonal resistance coefficient. The strong 
intrinsic spin-orbit coupling present in the host semicon- 
ductor valence band can account for the measured mag- 
nitude and sign of the AHE. 

There arc of course many other relevant experiments 
which arc not higlightcd above which must be taken into 
account. However the interpretation of some or most 
of them are still being debated, since they seem to de- 
pend critically on the annealing procedures and can be 
interpreted in several ways. The above results are con- 
sidered to be general and well established and will likely 
not change in the newly grown samples as the materials 
development progresses. 



II. THEORETICAL MODELS OF DMS 

The modeling of a collective behavior of interacting 
electrons is a complicated task and, in many instances, 
must be guided by experimental evidence of the low en- 
ergy degrees of freedom in order to obtain a correct 
minimal model which will capture the observed effects 
and will make useful predictions. With this in mind, 
there are typically three approaches, somewhat comple- 
mentary, used to describe DMS systems: first principles 
density-functional theory (DFT) and microscopic tight- 
binding models, effective Hamiltonian models, and lattice 
models. 

DFT is an important tool for studying microscopic ori- 
gins of ferromagnetism through calculations of electronic, 
magnetic, and structural ground-state properties. A 
local-density-approximation (LDA) of the DFT. com- 
bined with disorder-averaging coherent-potential approx- 
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imation (CPA) or supercell approach, has been used suc- 
cessfully to address physical parameters of (III,Mn)V 
DMS that are derived from total-energy calculations, 
such as the lattice constants,^'' and formation and bind- 
ing energies of various defects. ^^'^^ However, in Mn- 
doped DMS, LDA fails to account for strong correlations 
that suppress fluctuations in the number of electrons in 
the d-shell. As a result, the energy splitting between 
the occupied and empty d-statcs is underestimated which 
leads to an unrealistically large d-statc local DOS near 
the top of the valence band and to an overestimate of 
the strength of the sp-d hybridization. ^^'"^^ This short- 
comings have been corrected recently by LDA-I-U and 
self-interaction corrected LDA schemes, which account 
for correlations among Mn 3d electrons, and have been 
used to obtain more realistic energy spectra and show 
agreement with the experimental observation that the 
valence holes have mostly As 4p-character.'^^ 

A practical approach, that circumvents some of the 
complexities of this strongly-correlated many-body prob- 
lem is the microscopic tight-binding (TB) band-structure 
theory. Within the model, local changes of the crystal 
potential at Mn and other impurities are represented by 
shifted atomic levels. A proper parameterization of these 
shifts, of the Hubbard correlation potential that favors 
single occupancy of the localized d-orbitals, the Hund 
potential forcing the five d-orbital spins to align, and of 
the hopping amplitudes between neighboring atoms pro- 
vides correct band gap for the host IH-V semiconductor 
and an appropriate exchange splitting of the Mn d-lcvels. 
The TB model is a semi-phenomenological theory, how- 
ever, it shares the virtue with first principles approaches 
of treating disorder microscopically. The decoherence of 
Bloch quasiparticle states or effects of doping and dis- 
order on the strength of the sp-d exchange coupling and 
effective Mn-Mn interaction are among the problems that 
have been analyzed using this tool.'^'^"^^ 

In the metallic regime, where the largest critical fer- 
romagnetic temperatures are achieved ( for doping levels 
above 1.5%^°), semi-phenomenological models that are 
built on Bloch states for the band quasiparticles, rather 
than localized basis states appropriate for the localized 
regime,**^ provide the natural starting point for a model 
Hamiltonian which reproduces many of the observed ex- 
perimental effects. Recognizing that the length scales 
associated with holes in the DMS compoimds are still 
long enough, a k • p envelope function description of 
the semiconductor valence bands is appropriate. Since 
for many properties, e.g. anomalous Hall effect and 
magnetic anisotropy, it is necessary to incorporate in- 
trinsic spin-orbit coupling in a realistic way, the six- 
band (or multiple-band, in general) Kohn-Luttinger (KL) 
k • p Hamiltonian that includes the spin-orbit split-off 
bands is desirable. '''^^ The approximation of using the 
KL Hamiltonian to describe the free holes is based pri- 
marely in the shallow acceptor picture demonstrated by 
the experiments^^"^^ in (Ga,Mn)As and (In,Mn)As and 
must be re-examined for any other DMS materials that 



this model is applied to. 

Besides the KL Hamiltonian parameters of the host 
HI-V compound which have long been established,^^ the 
phenomenological part of the strategy follows from as- 
serting, rather than deriving, the localized nature of the 
Mn d-orbital moments and from parameterizing the sp- 
d hybridization by an effective exchange constant Jpd. 
The localization assumption is again verified by electron 
resonance expcrimcnts^''^^^^ and the value of Jpd is ob- 
tained from resistivity measurements in the paramag- 
netic regime^^ and MCD measurements.^^ Hence, the ef- 
fective Hamiltonian considered within this model is 

H = Hkl + Jpd S/ • s(r)5(r - Ri) + Hdis, (1) 

where H-kl is the six-band (multiple-band) Kohn- 
Luttinger (KL) k • p Hamiltonian,^'^ the second term is 
the short-range antiferromagnetic kinetic- exchange inter- 
action between local spin S/ at site R/ and the itinerant 
hole spin (a finite range can be incorporated in more real- 
istic models), and Tidis is the scalar scattering potential 
representing the difference between a valence band elec- 
tron on a host site and a valence band electron on a Mn 
site and the screened Coulomb interaction of the itinerant 
electrons with the ionized impurities. 

Several approximations can be used to vastly sim- 
plify the above model, namely, the virtual crystal ap- 
proximation (replacing the spatial dependence of the lo- 
cal Mn moments by a constant average) and the mean 
field theory description.^'^ In the metallic regime, the 
disorder can be treated by a Born approximation or 
by more sophisticated, exact-diagonalization or Monte- 
Carlo methods. The effective Hamiltonian in Eq. 1 
allows to use standard electron-gas theory tools to ac- 
count for hole-hole Coulomb interactions. This enve- 
lope function approximation model is also suitable for 
studying magnetic semiconductor heterostructures, like 
superlattices, quantum wells and digitally doped layers. 
The validity of such semi-phenomenological Hamiltonian, 
which does not contain any free parameters, must be 
confirmed ultimately by experiments. Its accurate de- 
scription of many thermodynamic and transport proper- 
ties of metallic (Ga,Mn)As samples, such as the mea- 
sured transition temperature, ^'^■^^ the anomalous Hall 
effect, ^^'^^ the anisotropic magnetoresistance,*^'^^, the 
magneto-crystalline anisotropy, the spin-stiffness"'''*, 
the ferromagnetic domain wall widths, the magnetic 
dynamic damping coefficients,^^ and the magneto-optical 
properties, ^'''^'^^ has proven the merit of this effective 
Hamiltonian approach. 

One has to keep in mind, however, that as any semi- 
phenomenological model it may fail to capture the cor- 
rect physics that leads to the ferromagnetic phase in 
some materials or in a certain range of parameter val- 
ues. Such models can only be verified by careful com- 
parison with experiments and tested through their pre- 
dictions and agreement with experimetnts. For exam- 
ple, Mn-doped nitride and phosphide compounds or in- 
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sulating DMS samples with low concentration of Mn ions 
require a theoretical description that goes beyond the 
picture of the host band quasiparticles that arc weakly 
hybridized with the localized Mn d-electrons. Particu- 
larly nitrates are not believed to be well modeled by this 
semi-phcnomenological Hamiltonian since Mn is a deep 
acceptors in this case and and charge fluctuations on the 
Mn d-levels may play important role. 

There has also been theoretical work on (lll,Mn)V 
DMS materials based on a still simpler model where 
holes are assumed to hop only between Mn acceptor 
sites, where they interact with the Mn moments via phe- 
nomenological exchange interactions. These models 
have the advantage of approaching the physics of the 
insulating dilute Mn limit, and can also be adapted to 
include the holes that are localized on other ionized de- 
fects besides the Mn acceptors through dynamical mean 
field (or CPA) techniques. The free-parameter nature of 
this phenomenological approach and their oversimplified 
eletronic structure allows to make only qualitative predic- 
tions, however, and the models are also not appropriate 
for studying the high Tc metallic samples. 



III. TRANSPORT PROPERTIES OF DMS 
SYSTEMS 

The different transport coefficients of DMS and their 
magnetic, temperature, and material composition depen- 
dance have been the most important and widely used 
characterization tools of DMS. Besides the diagonal re- 
sistivity which indicates metal-insulator transitions and 
possible critical behavior at the ferromagnetic transition 
temperature, other material transport properties such as 
anisotropic magnetoresistance, anomalous and ordinary 
Hall effect, and giant magnetoresistance have been used 
to both characterize and test the different theoretical 
models of DMS materials. In this section we consider the 
diagonal conductivity general features, how the exchange 
coupling between the localized moments and free carriers 
is extracted within the paramagnetic regime, the dop- 
ing and carrier concentration conductivity dependence 
at low temperatures, the anisotropic magnetoresistance, 
and the anomalous Hall effect. Throughout, we will fo- 
cus on the comparison of theoretical models (mainly the 
semi-phenomenological effective Hamiltonian model rele- 
vant to the metallic regime) to the different experimental 
observations. 



A. General features of resistivity in DMS 

DMS materials exhibit an insulating or metallic be- 
havior (defined by the resistivity in the limit of zero tem- 
perature) depending on its doping level and post-growth 
annealing procedures. In as-grown samples, metallic be- 
havior is typically observed for a range of 2-5% Mn dop- 
ing and an insulating behavior for higher and lower dop- 



ing than this range."'''' In addition to this metal-insulator 
quantum transition, the resistivity as a function of tem- 
perature typically exhibits a peak or shoulder near the 
ferromagnetic transition temperature for both insulating 
(peak) and metallic (shoulder) samples. ^^"^^ The non- 
monotonic behavior near Tc is typically associated with 
critical-scattering but so far no theory has been devel- 
oped which explains such behavior in a qualitative or 
quantitative way. Typically the Fisher-Langer theory of 
correlated fiuctuations is invoked, however, it predicts an 
infinite derivative of the resistivity at Tc, which is clearly 
not the case in any studied (Ga,Mn)As DMS sample. 
There exists also a drastic reduction of the resistivity 
upon annealing, associated with the increase of the car- 
rier concentration and to a lesser extent a reduction of 
the disorder scatterers.^^'^®"^^ The on-set of the metal- 
insulator transition at 1.5% Mn doping is close to the 
Mott insulator limit of a doped semiconductor similar to 
Si:P and the optimally annealed samples remain metallic 
throughout the whole range of Mn concentrations above 
1.5%.4o 

The number of research groups involved in the mate- 
rials growth process, each trying a slightly different an- 
nealing process, has increased over the past few years and 
with it a dramatic increase of carrier concentration and 
conductivity has taken place, as shown in Fig. 1. At 
the same time the Tc has also increased in accordance 
with the mean field theory prediction that Tc oc xp^^^, 
where x is the Mn concentration and p is the carrier con- 
centration. The shoulder in p observed in the most re- 
cent optimally annealed samples near Tc has been par- 
tially explained theoretically in terms of the variation of 
the Fermi surface and the transport scattering time as- 
sociated with the ferromagnetic to paramagnetic phase 
transition. ^"^ However, there is no model at present that 
reproduces fully the behavior observed in p as a function 
of temperature near Tc in the metallic regime for most 
samples. There has also been theoretical progress in un- 
derstanding the low temperature regime*^ (Sec. Ill C) 
and the role of scattering off magnetic impurities in the 
high temperature paramagnetic regime (Sec. IIIB).^^'^® 



B. Scattering off the kinetic-exchange potential in 
the paramagnetic regime 

The contribution from magnetic impurities to the re- 
sistivity behavior observed above Tc can be understood 
assuming scattering from the kinetic-exchange term in 
the Hamiltonian 

Wfe_e = Jpd^<5(r-R/)S/-s. (2) 
/ 

The corresponding contribution to the resistivity is given 

by32,33 

Ps = 27r2-^^^^i^J5^^iVMn[2x±(T,B) +X||(T,S)], (3) 



6 




X [Mn %] 



FIG. 1: Conductivity a measured at 4 K or lower vs. Mn 
percentage concentration for (Ga,Mn)As. The symbols corre- 
spond to: squares to FIX Ref. 24, circles to Ref. 64, triangles 
to Ref. 65, and dimonds to Ref. 59. 



where kp is the carrier Fermi wave vector, h is the Planck 
constant, Nmu the density of substitutional Mn, m* is 
the effective mass of the carrier, and e is the electron 
charge. X-L — MjB and xii = dM/dB are the trans- 
verse and longitudinal magnetic susceptibilities. Using 
transport data, the susceptibility can be determined from 
the Hall resistivity due to the dominant contribution 
from the anomalous Hall effect which is proportional to 
the magnetization.^^ The analysis of magnetoresistance 
data above Tc gives Jpd = 55± lOmeVnm'^. This result is 
in agreement with optical MCD data (see Sec. IV). We 
also note that the initial expression used to analize the 
magnetoresistance data which replaces the susceptibility 
factor by [S'(S'-|- 1) — (S)^], neglected the correlations be- 
tween neighboring Mn spins and overestimated Jpd by a 
factor of 3 in disaggrement then with the MCD measure- 
ment of Jpd-^^ 



C. Boltzman transport theory of the zero 
temperature conductivity in DMS 

The zero temperature conductivity of metallic DMS 
samples can be obtained from the effective Hamiltonian 
(Eq. 1) and by accounting for disorder scattering per- 
turbatively. The valence band holes interact with ran- 
domly located spins of substitutional Mn impurities via 
the kinetic-exchange interaction, and with randomly lo- 
cated ionized defects and each other via Coulomb inter- 
actions. At zeroth order, the interactions are replaced 
by their spatial averages, so that the Coulomb interac- 
tion vanishes and hole quasiparticles interact with a spa- 
tially constant kinetic-exchange field. The corresponding 



mean-field Hamiltonian for the itinerant holes reads 

Ha - Hkl + JpdNMnSn • s , (4) 

where fl is the orientation of fully polarized substitu- 
tional Mn local moments and s is the envelope-function 
hole spin operator.*^ Using the eigenstates of the Hamil- 
tonian in Eq. 4, the first order Born approximation 
of the elastic scattering rate, and the relaxation-time- 
approximation solution to the semiclassical Boltzmann 
equation, the diagonal dc conductivity tensor along one 
of the cube edges of the host lattice can be written as:'*'* 

where j is the quasiparticle elastic scattering rate, n 
and k are the band and wavevector indices, ^ are the 
eigenstates of the Hamiltonian (4), and Ep is the Fermi 
energy. 

The Born approximation estimate of the transport 
weighted scattering rate from substitutional Mn impu- 
rities is given by: 




X '5(£;„,,--i?„,s,)(l-cos0,-,,,), (6) 
with the scattering matrix element, 

<'5 = Jp<iS{zjn-Az^r,,) 



Here thost is the host semiconductor dielectric constant, 
|z^^) is the six-component (multi-component) envelope- 
function eigenspinor of the unperturbed Hamiltonian 
(4), and the Thomas-Fermi screening wavevector qtf = 
y/e^DOS{EF)/{ehosteo), where DOS{Ef) is the density 
of states at the Fermi energy.^** This model incorporates 
the fact that the transport properties of these materials 
are not solely determined by the scattering from substitu- 
tional Mn impurities and allows explicitly for scattering 
from compensating defects, which have been seen to play 
a key role in the resistivity through the post-growth an- 
nealing as discussed earlier. ^■^'^"^^ As-antisite defects are 
non-magnetic double-donors and contribute to scatter- 
ing through a screened charge Z = 2 Coulomb potential. 
The double-donor Mn interstitials^'*'^^ are unlikely to be 
magnetically ordered and can also be modeled as a Z = 2 
screened Coulomb potential.*^ 

Assuming a parabolic-band disperssion for majority 
heavy-hole states the kinetic-exchange scattering contri- 
bution to the scattering rate can be estimated by, Tpd = 
iNMn)JpdS'^rn*y/2m*EF/i4:nfi^). The Mn and As- 
antisite Coulomb scattering leads to scattering rate Tc 
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FIG. 2: Conductivity for measuring current and magnetiza- 
tion directed along the x-axis in the plane of the (Ga,Mn)As 
film as a function of the total hole density. These results were 
obtained for a GaAs semiconductor host doped with 6% Mn. 
For typical hole densities the current is carried mostly by the 
majority heavy-holes. 



FIG. 3: Experimental (filled symbols) AMR coefficients and 
theoretical data obtained assuming As-antisite compensation 
(open symbols) and Mn-interstitial compensation (semi-filled 
symbols) of AMRop and AM Rip. After Ref. 52 



given by the Brooks-Herring formula.''"' For (Ga,Mn)As, 
taking a heavy-hole effective mass m* = 0.5me, p = 
0.4 nm~^ and Mn doping x = 5%, these estimates give 
HTpd ~ 20 meV and HTc ~ 150 meV. A full numer- 
ical six-band calculations is consistent with these esti- 
mates, and predicts that the Coulomb contribution to 
the elastic scattering rate is several times larger than the 
kinetic-exchange contribution for typical chemical com- 
positions. Note that hese estimates give an immedi- 
ate check on the assumption of the theory itself, since 
even in the heavily doped and compensated (Ga,Mn)As 
DMS, the lifetime broadening of the quasiparticle (HT) 
is smaller than the valence band spin-orbit coupling 
strength {Ago = 341 meV) and the typical Fermi energy. 

Fig. 2 shows a XX, calculated numerically using the six- 
band Kohn-Luttinger model and Eqs. (5) and (6), for a 
fully strained Mno.06Gao.94As sample. The substrate - 
DMS lattice mismatch, eo = {asut — cidms)/(idms, is 
between -0.002 and -0.003 in this casc.2.65,66 rj,^^ j^^g^, 

lute conductivities predicted by this model are reason- 
ably consistent with experiment. ^^'^^'^^ The dissagree- 
ment for lower Mn concentrations {x < 4%) of the the- 
oretical conductivities is most likely due to some com- 
bination of inaccuracy in the scattering amplitude esti- 
mates, unaccounted sources of scattering, and, especially 
at small x, coherent scattering effects that eventually lead 
to localization observed as an upturn in p at the lowest 
temperatures. 

In addition to the AHE (Sec. HID), strong spin- 
orbit coupling in the semiconductor valence band leads 
also to anisotropies in the longitudinal transport coef- 
ficients. In particular, when the magnetization M is 
rotated by applying an external magnetic field stronger 
than the magneto-crystalline anisotropy field the in-plane 



conductivity changes. ■ Fig. 3 shows the theoretical 
and experimental AMR coefficients, AMRop ~ [o'xxi^ \\ 
z) - axx{^ II x)]/(yxx{y^ II x) and AMR,p = [axx{^ || 
y) - CTxx(M II i)]/crxx(M II x), for the seven (Ga,Mn)As 
samples. ^^'^^'^^ Here z is the growth direction. Results 
of the two disordered system models, one assuming As- 
antisite and the other one Mn-interstitial compensation, 
are plotted in Fig. 3. As in the AHE case, the theo- 
retical results are able to account semi-quantitativcly for 
the AMR effects in the (Ga,Mn)As DMS, with somewhat 
better agreement obtained for the model that assumes 
Mn-interstitial compensation, which confirms indirectly 
the experimental finding that the compensating defects 
are dominated by Mn-interstitials.^^ Note that although 
the magnitude of the conductivities tend to be overesti- 
mated, the magnetotransport effects are relatively insen- 
sitive to scattering strength, rcfiecting instead the strong 
spin-orbit coupling in the valence band of the host semi- 
conductor as compared to the Fermi energy.^^. 

A large anisotropic magnetoresistance effect is ob- 
served experimentally also in the off-diagonal symmetric 
conductivity component {axy = +(yyx) when magnetiza- 
tion direction changes from parallel to x- or y-axis, where 
Cxy = cFyx = 0, to a general orientation in the x-y plane 
tilted from the two transport measurement axes, where 
o'xy = <^yx 7^ 0.^* The effect is called sometimes a 'planar 
Hall effect' but we emphasize that the true Hall response 
requires the symmetry a^y = —(^yx- The anomalous con- 
tribution to this asymmetric off-diagonal transport coef- 
ficient is discussed in the following section. 
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D. Hall resistivity and the Anomlaous Hall effect 

One of the most important characterization tools of 
magnetic materials is the anomalous Hall effect^''. The 
Hall resistance RhuH = Pxy/d of a magnetic thin film 
is empirically observed to contain two distinct contri- 
butions. The first arises from the normal Hall effect 
contribution which is proportional to the applied mag- 
netic field, -ff, and the second, called the anomalous 
Hall contribution, is observed to be proportional to the 
magnetization^*^ : 

RHaii = RoH + RsM, (8) 

where M is the magnetization perpendicular to the thin 
film surface, and Rq and Rs are the ordinary and anoma- 
lous Hall coefficients respectively. The ordinary Hall ef- 
fect is due to the Lorenz force on the carriers and is used 
to measure the concentration and nature of the free car- 
riers since i?o is linearly proportional to the free carrier 
density and its sign determines whether they are elec- 
trons or holes. On the other hand, the anomalous Hall 
effect is a direct consequence of the presence of spin-orbit 
coupling in the system. In many instances, such as the 
DMS for example, the anomalous contribution is much 
greater than the ordinary Hall effect, making the carrier 
density hard to measure since large magnetic field must 
be used to extract the linearity in the field of the Hall 
resistance. However, the overwhelming of the Hall co- 
efficient hy Rs allows one to utilize the anomalous Hall 
effect as an indirect measure of the magnetization. Em- 
pirically Rs is observed to depend quadratically or lin- 
early on the diagonal resistivity, p, depending on the ma- 
terial. For example, for the case of Fe^° and GaMnAs''^ 
the dependance on the resistivity is quadratic. These 
two resistivity dependencies of Rs correspond to differ- 
ent possible origins of the anomalous Hall effect which we 
discuss below. However, in spite of its wide use, it is sur- 
prising how much confusion has followed over the past 
six decades. Although much progress has taken place 
over the past few years with the creation and study of 
new materials, no single theory has been able to system- 
atically and convincingly explain the constant puzzle of 
side-jump vs. intrinsic mechanisms (sec below) although 
it is at present a feasable but challenging undertaking. 

1. The embroiled history of the anomalous Hall effect 

Before discussing the theoretical description of the 
anomalous Hall effect in DMS, it is appropriate to make 
an excursion to the history of the field in order to un- 
derstand clearly the origin of the controversies and the 
strengths and weaknesses of the different arguments put 
forth by different researchers. Often such summaries are 
done in two or three sentences giving an oversimplified 
view of the controversy and the different points of view, 
some of which have been challenged in later literature 
but still remain in the general folklore. 



The first detailed theoretical consideration of the 
anomalous Hall effect was given by Karplus and 
Luttinger^^, where they considered the problem from a 
perturvative point of view (with respect to an applied 
electric field) and obtained a contribution to the Hall 
conductivity in systems with spin-orbit coupled Bloch 
states given by the expression (Eq. 2.16 and 2.17 in Ref. 
71) 




where np is the Fermi occupation number of the Bloch 
state \u^j:). This expression is obtained by ignoring im- 
purity scattering (i.e. clean limit) and can also be ob- 
tained from the Kubo linear response theory in this clean 
limit (see below). One of its immediate successes was 
the prediction that Rs is proportional to in agree- 
ment with many materials, Fe being the primary exam- 
ple. Immediately after Karplus and Luttinger's pioneer 
work, Smit^^ considered the problem of impurity scatter- 
ing within a Boltzman formalism using a model Hamilto- 
nian of an electron gas, i.e. plain waves, without intrinsic 
spin-orbit coupling and a scattering potential containing 
both direct and spin-orbit coupled terms: 

^=£ + ^(^)+2^^^^^^ 

where the impurity potential is V{r) = for r > R 
and V{r) = Vq for r < R. Within the Boltzman for- 
malism, he obtained a second type of contribution to 
the Hall conductivity due to the asymmetric scattering 
from impurities originating from the spin-orbit coupling 
term. This so called skew-scattering mechanism predicts, 
in contrast to the intrinsic effect found by Karplus and 
Luttinger, that Rs oc p (i.e. density of scatterers) and 
dependens on the type and range of the scattering poten- 
tial. It was also noted that the magnitude of the effect 
within such simplified model was many orders of magni- 
tude lower than the observed magnitude of the AHE but, 
as it was shown later, the origin of this deficiency is the 
lack of inrtinsic spin-orbit coupling present in the crys- 
tal Bloch states. However, Smit also critiziced Karplus 
and Luttinger's result claiming the intrinsic contribution 
vanishes. Such claim (see below) is now understood to 
be unjustified but has remained lingering in the AHE 
literature, generating much confusion. 

Karplus and Luttinger''^ did not consider impurity 
scattering, hence omitted the skew-scattering contribu- 
tion in their calculation. To remedy this, Kohn and 
Luttinger'''^ developed a consistent (and somewhat coum- 
bersome) treatment to obtain the transport coefficients 
in the presence of spin-orbit coupling and disorder scat- 
tering based on a density matrix expansion. Such ap- 
proach is equivalent to the Kubo linear response formal- 
ism but technically more difficult when applied to real 
physical systems. Within this formalism the expansion 
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is done to first order in spin-orbit coupling and impurity 
scattering strengths. They obtained, besides the previ- 
ous Karplus and Luttingcr result (Eq. 135 in Ref. 73, 
however note that it has the wrong sign later corrected 
in Ref. 74), formal expressions for additional contribu- 
tions due to impurity scattering. Applying this formal- 
ism to a model with uncorrelated impurity potential and 
small scattering strength, Luttinger'^^ obtained, besides 
the Smit's skew-scattering term, the same final expres- 
sion as Karplus and Luttinger but with the opposite sign 
(Eq. 3.26 in Ref. 74). Somehow, in addition to the 
intrinsic Karplus and Luttingcr term, there seemed to 
be an additional contribution from scattering in the high 
mobility limit giving a contribution identical to the in- 
trinsic one but with a factor of —2 in front of it, hence 
the flip of sign of the final expression. No explanation 
was given of why such a scattering contribution does not 
depend in any way on the scattering potential but sim- 
ply on the electronic structure of the system. However, 
Luttingcr^'* showed inequivaqucly that the cancelation 
argued by Smit does not take place. In a slightly different 
formalism Adams and Blouht^^ agreed with Luttinger's 
results and pointed out that Smit's error occurs from an 
inconsistency involving a change of representation of his 
density matrix when calculating the current. 

A few years later Bcrgcr^*' made his first contribution 
to the problem introducing the idea (contained withing 
Luttinger's formalism) of the side-jump scattering mech- 
anism. He considered the problem of wavc-pakt^t scatter- 
ing off an impurity potential as in Eq. 10, showing that in 
addition to skew scattering there exists a side-step type of 
scattering. This so called side-jump, as in the Luttinger 
formalism, was independent of the scattering potential, 
leaving one with the strange notion that an exterinsic 
scattering mechanism resulted in an intrinsic property. 
After putting this contribution in a Boltzman type of 
formalism he obtained a result which was four orders of 
magnitde lower than the Karplus and Luttinger intrinsic 
result. He then argued that the intrinsic spin-orbit coulp- 
ing of the system increases this side-jump scattering by 
such factor. In his case the sign of the effect is not di- 
rectly considered and his result, as already mentioned, is 
independent of the scattering potential itself. A similar 
type of formalism was put forth by Lyo and Holstein''^ 
claiming similar results but ignoring, for some unknown 
reason, intrinsic non- vanishing contributions of Karplus 
and Luttinger. 

It is then that Leroux-Hugon and Ghazali^* pointed 
out the resolution of the missing order of magnitude 
of the skew-scattering. By looking simply at scattering 
of Bloch electrons from ionized impurities Vton {r) alone 
and going to second order in the Born approximation 
within the collision term of the Boltzman transport the- 
ory, they obtained a skew-scattering contribution with a 
much higher magnitude than the one coming from the 
spin-orbit coupling of the impurities themselves. Hence, 
it became evident that it is the presence of intrinsie spin- 
orbit coupling through the crystal potential which is ul- 



timately responsible for the magnitude of all the anoma- 
lous Hall effect contributions observed and a system with 
weak or no intrinsic spin-orbit coupling will not exhibit 
a detectable contribution to the RhuU ■ 

To elucidate this confusing and farraginous atmo- 
sphere, Nozieres and Lewiner^^, following the formalism 
introduced by Fivaz^° in terms of effective Hamiltoni- 
ans of semiconductors, developed a theory applicable to 
narrow gap semiconductors equivalent to Luttinger's^^ 
but with the final result in being more physically trans- 
parent. Going to linear order in impurity scattering and 
spin-orbit coupling and combining them with a Boltzamn 
type of approach, they obtained different contributions to 
the Hall current and different behavior when considering 
opposite orders of limits of w ^ and r ^ oo, where to 
is the frequency of the applied electric field and r is the 
quasiparticle lifetime (see Table I and Eqs. 59, 60a, and 
60b in Ref. 79). Besides the skew scattering contribu- 
tion which depends on the type of scattering potential, 
in the limit lot » 1 (r ^ oo first), considered by some 
authors to be the weak scattering limit (see Refs. 81 
and 82), the result is that of Karplus and Luttinger.''^ 
On the other hand, for lot << 1 (cj ^ first), the oppo- 
site sign expression of Luttinger^^ is recovered. One must 
note however that this result seems to be very specific to 
the model (single band) and to the meaning of the limit 
T ^ 00 as admitted by the authors themselves. '^^■''^ 

After this Berger and Smit exchanged salvos in a se- 
ries of confusing comments*^"^^ in which Smit's old argu- 
ment, that there is no such a thing as side-jump or intrin- 
sic contribution, seems to have been finally put to rest. It 
was after this that Chazalviel*^ , using Nozieres-Lewier- 
Luttinger^^'™ and Leroux-Hugon and Ghazalis^^ results, 
attempted to compare phenomenologically the AHE the- 
ory to n-InSb and n-Ge experimental data with some 
success. Within such models, the relative importance of 
the sidc-jump/intrinsic contribution c^^*"* (which are 
the same magnitude but of opposite sign depending on 
the limits^^'*^'^^) to the skcw-scattring contribution a^Jjj 
depends on the nature of the scatterers (only through the 
skew-scattering dependence). For scattering from ionized 
impurities one has''^'^^ 



sj /int 



cN 



(11) 



with N/p being the ratio of the total number of ion- 
ized impurities and carrier density, is the average 
distance between carriers in units of Bohr radious, I is 
the mean free path, and c ^ 10, varying slightly with 
scattering length. For short range scattering potential 
considered by Luttinger''^ and Nozieres and Lewiner^^, 
V{r) = Vo6{f-r-l): 



sj /int 



'aH 



Tl\Vo\D{EF)kFl 



(12) 



where D{Ep) is the density of states at the Fermi energy 
and kp is the Fermi wave- vector. These estimates can be 
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a useful first guess at which mechanism dominates in dif- 
ferent materials but one must keep in mind the simplicity 
of the models used to estimate such ratio. 

In spite of the theory being incomplete, and in some 
cases self-contradictory and puzzling, the theoretical de- 
velopments remained somewhat dormit until recent inter- 
est in new ferromagnetic and spintronic materials. These 
emerging field have demonstrated the need for a better 
theoretical understanding of the AHE which has con- 
stantly been used to charactcizc ferromagnetic materials. 

We emphasize that there is, at this stage, no con- 
troversy regarding the skew-scattering contribution. In 
most materials such contribution, although present, is 
minor and in most ferromagnetic materials the quadrac- 
tic dependence of Rs is observed even at the lowest tem- 
peratures. However, if somehow, one could artificially 
turn down p one would eventually reach a regime where 
such contribution is dominant, a fact which very few 
people dispute since the origin of such contribution is 
rather transparent and reasuringly extrinsic in nature. 
However, the flip of sign of the expressions derived by 
Luttinger'''* and Nozieres and Leweir*"^ within the sim- 
ple one band model in the different order of limits of 
LOT, remains quite disturbing and a source of continuous 
debates. After all, these results are only obtained in a 
limit where, siipposcdly, the skew scattering would al- 
ways dominate. As we will see below, it seems that in 
many materials experimental comparisons, the ultimate 
test of validity of any theory, are more consistent with 
the intrinsic Karplus and Luttinger ''^ contribution than 
the reversed sign one of Luttinger''*; the sign of the effect 
being a simple thing to check. 

A possible resolution to the problem would be to per- 
form the Hall conductivity calculation within the Kubo 
linear response formalism, treating disorder and spin- 
orbit coupling on an equal footing. This can be achieved 
by starting from the Dirac representation and taking the 
weak relativistic limit after treating disorder in the usual 
fashion. In this way, any spin-orbit coupling contribu- 
tion comes naturally from the Dirac representation and 
therefore is automatically taken into account. On the 
other hand, if one starts from the Pauli Hamiltonian, an 
inmediate question arises: which vertex corrections cap- 
ture the different contributions to the anomalous Hall ef- 
fect? Such an approach was taken recently by Crepieux, 
Dugaev, and Bruno. ^''"^^ The work considers, however, 
a free electron gas model rather than a crystaline cnvi- 
roment, hence ignoring the intrinsic spin-orbit coupling 
effects. In spite of this simplification, the authors ob- 
tained many useful findings. They where able to pin- 
point which diagrams correspond to the skew and the 
side-jump scattering in the Pauli Hamiltonian by iden- 
tifying them from the corresponding ones arising natu- 
rally from the Dirac formalism. Notably, these diagrams 
within the Pauli Hamiltonian, are not the standard ones 
that one would take into account naturally. Their main 
result, Eq. 50 and 51 in Ref. 87,89, is in reasonable 
agreement with Berger's simpler treatment of this free 



electron gas problem. 

Taking into consideration the confusion generated by 
the possible relevance of the side-jump scattering mecha- 
nism, several researchers have chosen to focus instead on 
the original intrinsic contribution to the AHE proposed 
by Karplus and Luttinger''^ and ignore impurity scat- 
tering all together or simply include its effects through 
the Born approximation which introduces a finite quasi- 
particle lifetime. In DMS, e.g., this approach is partly 
motivated by the strong intrinsic spin-orbit coupling in 
the host valence bands that makes much of the above 
theoretical discussions inapplicable since they relied on 
the perturbation treatement of the spin-orbit coupling. 
One of the main attraction of the intrinsic AHE theory 
is the ability to do calculations in models with realis- 
tic electronic band-structures. This approach has been 
used, e.g., to analyze the AHE in layered 2D ferromag- 
nets such as SrRuOs, in pyrochlore ferromagnets^""^'^ , 
in the collosal magnetoresistance manganites,^'' in Fe,^^ 
in DMS,^^ and as a natural starting point to address 
infrared magneto-optical effects such as the Kerr and 
Faraday effects^'' (see Sec. V). The application of this 
approach to these different materials and their succesful 
comparison to experiments^^'^^'^^'^°~^^'^^'^^ is perhaps 
the most pogniant criticism to the old theories regard- 
ing the side-jump scattering as fundamental. 



2. Anomalous Hall effect in DMS 

The anomalous Hall effect in DMS has been one of 

the most fundamental characterization tools since it is 
the simples way of detecting the ferromagnetic state of 
the system at a given temperature. The original dis- 
covery of ferromagnetism in (III,Mn)V seminconductors 
was established by measuring the AHE in both the high 
and low temperature regimes.^ The comparison with 
remenant magnetization measurements using a SQUID 
magnetometer confirmed, e.g., that AHE measurement 
can be used to determine ferromagnetic critical temper- 
ature.^ 

Recently Jungwirth, Niu, and MacDonald^^ reintro- 
duced the original Karplus and Luttinger Eq. 9, point- 
ing out that the intrinsic contribution to the AHE is pro- 
portional to the Berry phase acquired by a quasiparticle 
wave function upon traversing closed path on the spin- 
split Fermi surface. They applied the theory to metallic 
(III,Mn)V materials using both the 4-band and 6-band 
k-p description of the valence band electronic structure"^^ 
and obtained results in a quantitative agreement with the 
experimental data in (Ga,Mn)As and (In,Mn)As DMS. In 
DMS systems, the estimate given in Eq. 11 gives a ratio 
of intrinsic to skew scattering contribution of the order 
of 50, hence the intrinsic contribution in these systems 
is likely to dominate. Consistently, Edmonds et al.^^ 
found that in metallic DMS systems Rg oc p^. In a follow 
up work^^ a more carehil comparison of the theory, refor- 
mulated within the Kubo formalism, was done in order to 
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account for finite quasiparticlc lifetime effects important 
only for quantitative but not qualitative comparison with 
the experimental data. Within the Kubo formalism the 
dc Hall conductivity for non-interacting quasiparticles is 
given by 



-■xy 



ie^h f dk -sr-^ 
J (2^ 



E E ,t 



{nk \px I n'k) {n' k\py\ nk) 
E ,t - E t + ihV 



(13) 



Looking at the real part of the dc Hall conductivity in 
the clean limit {hT 0), the delta- function contribution 
from the denominator vanishes due to the Fermi factor 
differences and we obtain 



-"xy 
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(14) 



Realizing that the dipole matrix elements con- 
sidered above can be written as {n'k\pa\nk) = 
{m./h){n'k\dH{k)/da\nk), Eq. 14 can be shown to be 
equivalent to Eq. 9. However, within the Kubo formal- 
ism, it is straightforward to account for the finite lifetime 
broadening of the quasiparticles within the simple Born 
approximation by allowing T above to be finite. The ef- 
fective lifetime for transitions between bands n and n', 
T„_„' = l/r„^„', can be calculated by averaging quasi- 
particlc scattering rates calculated from Fermi's golden 
rule including both screened Coulomb and exchange po- 
tentials of randomly distributed substitutional Mn and 
compensating defects as done in the dc Boltzman trans- 
port studies''^'^^ . In Fig. 4 the compensation is assumed 
to be due entirely to As-antisite defects. The valence 
band hole eigenenergies E^^ and eigenvectors \nk) in 
Eqs. (13)-(14) are obtained by solving the six-band Kohn- 
Luttinger Hamiltonian in the presence of the exchange 

field, h = NMnSJpdZ^'^ ■ Here Nmu = 4:x/<1dms 
Mn density in the Mn2;Gai__a;As cpilayer with a lattice 
constant udms, the local Mn spin S = 5/2, and the ex- 
change coupling constant Jpd = 55 meV nm~^. 

Fig. 4 demonstrates that whether or not disorder is 
included, the theoretical anomalous Hall conductivities 
are of order 10 cm"^ in the (Ga,Mn)As DMS with 
typical hole densities, p ~ 0.5 nm~^, and Mn concentra- 
tions of several per cent. On a quantitative level, disorder 
tends to enhance a ah at low Mn doping and suppresses 
AHE at high Mn concentrations where the quasiparticlc 
broadening due to disorder becomes comparable to the 
strength of the exchange field. The inset in Fig. 4 also 
indicates that the magnitude of the AHE in both models 
is sensitive not only to hole and Mn densities but also to 
the lattice-matching strains between substrate and the 
magnetic layer, cq = {asubstrate - aDMs)/aDMS- 




^ w w disorder, p=0.4 nm" 
e„=0 , 
" BQ clean, p=0.6 nm 

disorder, p=0.6 i^m"^ 



FIG. 4: Theoretical anomalous Hall conductivity of 
Mna:Gai_a;As DMS calculated in the clean limit (open sym- 
bols) and accounting for the random distribution of Mn and 
As-antisite impuriries (filled symbols), after Ref. 52. 



A systematic comparison between theoretical and ex- 
perimental AHE data is shown in Fig. 5.^^ The results 
are plotted vs. nominal Mn concentration x while other 
parameters of the seven samples studied are listed in the 
figure legend. The measured aAH values are indicated by 
filled squares; triangles are theoretical results obtained in 
the clean limit or for a disordered system assuming either 
the As-antisite or Mn- interstitial compensation scenario. 
In general, when disorder is accounted for, the theory is 
in a good agreement with experimental data over the 
fidl range of studied Mn densities from x = 1.5% to 
X = 8%. The effect of disorder, especially when assum- 
ing Mn-interstitial compensation, is particularly strong 
in the x = 8% sample shifting the theoretical (tah much 
closer to experiment compared to the clean limit theory. 
The remaining quantitative discrepancies between theory 
and experiment have been attributed to the resolution in 
measuring experimental hole and Mn densities. 

We conclude this section by mentioning other stud- 
ies of the anomalous Hall effect in regimes and geome- 
tries which we have not focused on here. A study within 
the non-metallic regime or hopping transport regime, not 
covered in any detailed within this review, has been done 
by Burkov et al.^"^ following the ideas utilized in the met- 
talic Berry's phase approach but generalized to hopping 
transport regime. Their results are in good agreement 
with experimental data in these regime which encom- 
pass samples with less than 1% doping. At tempera- 
tures above T^, a theory based on semiclassical Boltzman 
transport has suggested the presence of anomalous Hall 
voltage noise within this regime which could be a further 
test of the intrinsic AHE model. ^® 

A recent theoretical study by Bruno et al.^^ has sug- 
gested an alternative experimental system to test the 
topological Hall effect arising from the Berry phase sim- 
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A. visible and infrared absorption 
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FIG. 5: Comparison between experimental and theoretical 
anomalous Hall conductivities, after Ref 52. 



ilar to Eq. 9. Here the difference is that the system has 
a spatially -varying magnetization. Upon a local gauge 
transformation which makes the quantization axis ori- 
ented along the local magnetization direction, the prob- 
lem becomes that of a gauge vector potential correspond- 
ing to monopole fields. Such a system can be created 
by sandwiching a thin film (II,Mn)VI semiconductor and 
a lattice of iron nano-rods oriented perpendicular to the 
doped semiconductor epilayer. Future experiments might 
show whether such a device can be grown efficiently and 
the theory tested. 



IV. OPTICAL PROPERTIES OF DMS 

Optical properties are among the key probes into the 
electronic structure of materials and, in the case of semi- 
conductors, can also be utilized to change its proper- 
tics by photo-carrier generation and other means. Given 
the wide energy range that optical probes can attain, 
many different physics can be addressed by exploring 
the full spectrum of phenomena including the visible 
and infrared absorption, magneto-optical effects 
such as the magnetic circular dichroism,'^^'^''^"^''^ Ram- 
man scattering, photoemission,^^'^^'^^^'^^'^ and cy- 
clotron resonance. ^^'^"^^^ We focus in this review on the 
first two items of this long list. The motivation for the 
study of optical properties of materials is two fold: first, 
it gives a consistancy check on the different theoreti- 
cal models and can be used to measure directly phc- 
nomenological parameters such as the exchange coupling 
Jpd,^^"^ '^^^ and second, it can be utilized in applications 
such as magneto-optical memories, optical isolators, and 
circulators for optical communications. ^^'^ 



The simplest and most direct of the optical effects is 
light absorption, which probes the electronic structure 
through electron excitations between different bands or 
impurity states. An important measurement in semicon- 
ductors is the band edge absorption, usually obtained 
in the visible frequency range. This visible regime has 
been exploited more often using absorption from po- 
larized light in order to measure magneto-optical ef- 
fects such as the MOD and we will postpone its disuc- 
ssion to Sec. IV B. However, in the infrared regime, 
many of the meV physics is revieled and several re- 
markable experimental features have been observed in re- 
cent expcrments.'^-''^""'^"^'^"''"-^"''' These experiments, per- 
formed in thin film geometries, exhibit several common 
phenomena: (i) a non-Drude behavior in which the con- 
ductivity increases with increasing frequency in the inter- 
val between meV and 220 meV, (ii) a broad absorption 
peak near 220 — 260 meV that becomes stronger as the 
samples are cooled, and (iii) a broad featureless absorp- 
tion between the peak energy and the effective band gap 
energy which tends to increase at the higher frequencies. 

In order to understand these features, two theoreti- 
cal approaches have been employed mostly. The first is 
the scmi-phcnomcnological effective Hamiltonian model 
that focuses on the multi-band nature and on the in- 
trinsic spin-orbit coupling present in the host semicon- 
cutor in order to quantitatively understand these ma- 
terials and is able to capture many of the observed in- 
frared properties. The second approach, used in 
the context of lattice models, emphasizes the connec- 
tion with the localization physics present in the low dop- 
ing regime and simplifies the electronic structure by as- 
suming a single-band, seeking a more qualitative rather 
than quantitative understanding of these materials in the 
regime of interest. Both models arc important to attain 
a full understanding of DMS optical properties and we 
review the results of each in turn. 



1. Multiband effective Hamiltonian approach 

Within the effective Hamiltonian model (Eq. 1), the 
220 meV peak can be attributed to inter-valence-band 
transitions,^^ rather than to transitions between the 
semiconductor valence band states to a Mn induced im- 
purity band, although a combination of these contribu- 
tions will always be present. In thin film absorption mea- 
surements, for infrared wavelengths much larger than the 
width of the film, the real part of the conductivity is re- 
lated to the absorption coefficient by 



a{(j) = 2 



Re[a{uj)] 
Y + Yc 







(15) 



where Y and Yq are the admittances of the substrate and 

free space, respectively. For shorter wavelengths multi- 
ple refiections within the film must be taken into account 
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to obtain an cstiamte of the absorption coefficient. The 
conductivity tensor aa^isij-o) at T = can be evaluated us- 
ing the standard Kubo formula for non-interacting quasi- 
particles. Disorder is taken into account within the Born 
approximation by including the lifetime broadening of 
quasiparticle spectral functions in evaluating the Kubo 
formula. The effective lifetime for transitions between 
bands n and n', t„^„/, is calculated by averaging quasi- 
particle scattering rates calculated from Fermi's golden 
rule including both screened Coulomb and exchange in- 
teractions as in Sec. Ill C. Fig. 6 shows the ac conduc- 
tivities calculated for Ga.95Mn.05As at a series of carrier 
densities. 

An imporant concept in optical absorption spectra is 
the f-sum rule: 




'0 100 200 300 400 500 600 700 8o8 

CO [meV] 



F= / rfa;Re[a^^H] = — ("|-^^2q^I")- 

■'O ^ X 

(16) 

In this equation is a quasiparticle Fermi factor and 
d^Hxh/^^dk^ is the xx component of the k.p model in- 
verse effective mass operator. This sum rule is completely 
independent of the weak-scattering approximations but 
it is necessary to choose an upper cut-off for the fre- 
quency integral which creates a small uncertainity in the 
optical effective mass, defined by _F = ne^'pj^m.o.pi. As 
illustrated in the inset of Fig 6, for the case of cut-off 
frequency fuVmax = SOOmeV, the f-sum rule values of F 
evaluated from our weak scattering theory arc accurately 
linear in p over the entire range of relevant carrier den- 
sities. Disorder does have a small but measurable effect 
on niopt as illustrated in Fig. 6. The optical masses for 
GaAs, InAs, and GaSb DMS fcrromagnets with a 800, 
400, and 700 meV cutoffs are 0.25-0.29 nie, 0.40-0.43 
me, and 0.21-0.23 me respectively, the extremes of these 
ranges corresponding to the clean (lower) and disordered 
(upper) limits of our model. 

Free carrier concentration is difficult to determine ac- 
curately in DMS systems because the anomalous Hall 
contribution overwhelms the oridnary Hall coefficient 
and, in many instances, a high magnetic field is required. 
Hence, the f-sum rule can be used as an alternate tool 
to measure the free carrier concentration in metallic sys- 
tems where the assumed electronic structure may be a 
good approximation so niopt is given correctly by the es- 
timates above. On the other hand, accurate estimates of 
the carrier concentration p can be used to test the the- 
ory as well. Recent experiments^^'''^^''- in as-grwon sam- 
ples have been interpreted in both ways without clear 
distinction between the effective Hamiltonian picture or 
the impurity band picture since no measurements of the 
carrier concentration where available. Further measure- 
ments in the most metallic samples will serve as test of 
the simplifying assumptions within this model. 

One prominent feature that is at odds with experimen- 
tal data from the above model calculations is the relative 
magnitude of the lu ^ conductivity and the 220 meV 
conductivity peak. The reason is the impossibility of the 



FIG. 6: Optical conductivity Re[a(aj)] and absorption coef- 
ficient a(tj) for carrier densities from p = 0.2 to O.Snm^'* 
in tlie direction indicated by tlic arrow, for Gao.95Mno.05As. 
The dot-dashed line is an experimental absorption curve for a 
sample with a Mn concentration of 4% obtained from Ref. Hi- 
rakawa:2002. The inset shows the spectral weight F{p) eval- 
uated with a frequency cut-off of 800 meV including disorder 
(dashed) and in the clean Umit (solid) . After Ref. 45 

model to account for weak and strong localization effects 
from multiple-scattering effects invariably present in the 
DMS materials and which suppress te low-frequency con- 
ductivity. Exact diagonalization calculations still within 
the effective Hamiltonian model in a finite system size but 
treating the disorder effects exactly can fix this shortcom- 
ing as shown in Fig. 7 .^^ In this calculation the f-sum 
rule value differs from its Born approximation value by 
less than 10% for typical metallic carrier densities, hence 
justifying the use of niopt as a measuring tool of the car- 
rier concentration p. 

2. Lattice models and dynamical mean field theory approach 

Within the lattice models, dynamical mean-field- 
theory studies, ^^^'^^^ for a single-band model that neglect 
the spin-orbit coupling and the heavy-light degeneracy 
of a III-V semiconductor valence band, have sown that 
non-Drude impurity-band related peaks in the frequency- 
dependent conductivity occur in DMS ferromagnet mod- 
els when the strength of exchange interaction coupling 
is comparable to the band width. The starting Hamil- 
tonian of these theories is the generalized Kondo lattice 
model 

H = Hhost - ^ JSi- tljl{Ri)aal3tpl3{Ri) 

+ Wiji{Ri)ilJc.{Ri), (17) 

where Hhost describes carrier propagation in the host 

disordered semiconductor, here approximated by a sin- 
gle band Hamiltonian with semicircular density of states, 
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FIG. 7: Absorption and conductivity of a metallic sample 
computed using a Luttinger Kohn Hamiltonian with disorder. 
Here p = 0.33 nm"^, umu = 1 nm"* {x « 4.5%). After Ref. 
46 



the main band and the impurity band that forms when 
the carrier-Mn coupling becomes comparable or stronger 
to the main band width and spin-splitting. The high 
value of the exchange coupling required for this physics 
to apply is, however, not consistent with experimental 
Jpd values inferred from the visible-range MCD or mag- 
netotransport measurements. 

Monte Carlo simultations in closely related lattice 
models exhibit a similar conductivity behavior arising 
again from the presence of the impurity band.^^^ These 
theories point out the qualitative features present in the 
theoretical models as a function of the material parame- 
ters which, although fixed in most DMS materials, may 
be tuned by chemical engineering. It would be of inter- 
est to extend the calculation to a multi-band Hamilto- 
nian model which may lead to a more complex descrip- 
tion of the intermediate frequency regime. In addition, 
this extension would allow for a calculation of magneto- 
optical effects within this models which are not available 
at present. 



the second term describes the exchange coupling of the 
carriers to an array of Mn moments Si at positions Ri, 
and the third term is the scalar part of the carrier-Mn 
potential. At zero temperature, when all Mn moments 
Si are aligned, the carriers with spin parallel to Si feel 
a potential —J + W on each magnetic impurity site and 
anti-parallel carriers feel a potential J + W. Within this 
model, the disorder is treated in the dynamical mean 
field approximation, or equivalently in these systems the 
dynamical coherent potential approximation. The real 
part of the conductivity is given by 



a{n,T) = 



d^p f pcos9 



{2n) 



dw [/(w) - f{uj + Q)] 



n 



X ImG(p, w)ImG(p, u + 



(18) 



Typical results for the conductivity obtained within 
this model are shown in Fig. 3 of Ref. 118. Although 
the curves are similar to the ones obtained using the ef- 
fective Hamiltonian multi-band models, the origin of the 
peak at intermediate frequencies is quite different. In the 
lattice model the peak corresponds to transitions between 



B. Magneto-optical effects 

Magneto-optical effects, such as magnetic circular 
dichroism (MCD), Kerr effect, and Faraday effect, give 
further insight in ferromagnetic materials and add an adi- 
tional insight for the modeling of the electronic struc- 
ture. Absorption and reflection measurements in the 
visible range have been used to establish phenomeno- 
logical estimates for the p-d and s-d exchange cou- 
pling constants in DMS materials and are found to 
be in agreement with magnetic susceptibility measure- 
ments obtained from magneto-transport measurements 
in the paramagnetic regime. Measurements of 
magneto-optical coefficients on band energy scales pro- 
vide very detailed information about the influence of bro- 
ken time-reversal symmetry on itinerant electron quasi- 
particle states. For DMS materials the corresponding en- 
ergy scale for the heavily p-doped (III,Mn)V ferromag- 
nets is in the infrared. Within the effective Hamilto- 
nian model, the anomalous Hall effect theory^^ which 
has proven succesful when comparing closely with exper- 
imental results, can be easily extended to finite frequency 
range:^'' 



(f - - f -) 



e^h dk 

Im[(n'fc|p.|nfc)(nfc|p^|n'fc)](r2 „, +u,(E„^ - E^^p) - (E„g 



{{uj - E^^ + E^,^y + {hTr,,^.y){{E^^ - E^,^Y + {hT^,n'?) 

I 



(19) 



where r„^„' = (r„ -|- r„/)/2 and r„ are the golden rule lllC and IV A 1. Typical values for this anomlaous ac- 
scattering rates averaged over band n as shown in Sec. 
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FIG. 8: Anomalous ac-Hall conductivity crxy{i^) for x = 6% 
Mn concentration and p = 0.4 and 0.2nm~^, for spherical and 
non-spherical (band-warping) models. After Ref. 57 



Hall conductivity are shown in Fig. 8 

The Hall conductivity a^y must be non-zero in order 
to have non-zero magneto-optical effects, such as the 
Faraday and Kerr effect, but the full conductivity ten- 
sor must be known in order to connect the calculations 
with the measurable quantities. The calculations and 
measurements mentioned here are for the Faraday geom- 
etry (magnetic field parallel to light propagation axis) for 
transmission and the polar Kerr geometry (near normal 
incidence with the magnetic field along the light prop- 
agation direction) for reflection. The complex Faraday 
angle includes both real and imaginary terms. Op is 
defined as 



tan 9t. 



= OF 



(20) 



where t^x and txy are the complex transmission ampli- 
tudes for linearly polarized light and t± arc the total com- 
plex transmission amplitudes (with multiple scattering 
taken into account) for right and left circular polarized 
light. In typical measurements, 9f is small and therefore 
the small angle approximation tanflp ~ can be used. 
The real term (Re[0F] or Op) corresponds to a simple 
geometric rotation of the polarization vector about the 
direction of propagation. The imaginary term (Im[0F] or 
r]-p) relates directly to the ellipticity of the polarization. If 
the sample is axially symmetric along the magnetic field 
B, the transmittance tensor is diagonal when represented 
in the circular polarization basis. Therefore, changes in 
the incident polarization only depend on: (1) the rela- 
tive difference in the phase of left versus right circularly 
polarized light due to Re[0F], which leads to a rotation 
(circular birefringence or Faraday rotation, FR) in the 
linearly polarized incident fight; and (2) the relative dif- 
ference in the transmission of left versus right circularly 



polarized light due to Im[0F], which introduces ellipticity 
(circular dichroism, CD) to the linearly polarized inci- 
dent light. 

In the thin film approximation, the relationship be- 
tween Of and and the conductivity is given by: 



tan^F = 



_ txy 



^ (21) 



where Z ~ Zod/{ji+l), Zq is the impedance of free space, 
n is the substrate index of refraction, and d is the film 
thickness. In bulk, the angle of rotation per unit length 
traversed is 



An 



(22) 



where c is the speed of light and n is the index of refrac- 
tion of the substrate, in this case GaAs with n = y/lQ.9. 

CD or MCD (magnetic circular dichroism) is related 
to the difference between the optical absorption of right 
and left circularly polarized light, and is given by 



MCD = 



Re[CTa;j;(w)] 



lm[0F 



(23) 



where a± are the absorbances for left and right circularly 
polarized light. At visible range frequencies multiscatter- 
ing reflections must be taken into account since in that 
regime the wavelength is comparable to the typical epi- 
layer thickness. 

Whereas the complex Faraday angle describes the po- 
larization of a transmitted beam, the complex Kerr angle 
6k is applied in the same way for reflected radiation, with 



tan( 



(24) 



where r^x and r^y are the complex reflection amplitudes 
for linearly polarized light and r± are the total complex 
reflection amplitudes (with multiple scattering taken into 
account) for right and left circular polarized light. 

Fig. 9 shows the different magneto-optic effects pre- 
dicted for a concentration of x = 6% and p = 0.4 
nm~^. The Faraday rotation in this case is larger than 
the giant Faraday rotation observed in the paramagnetic 
(II,Mn)VI's at optical frequencies^^^ and should be read- 
ily observable in the highly metallic samples. 

Though magneto-optical measurements on II- VI DMS 
materials have been extensive, see for example Ref. 121, 
measurements on III-V DMS have been more limited, 
with almost no Faraday/Kerr measurements reported 
on these compounds in the mid-infrared (MIR, 100-500 
meV) range. The MIR range provides a critical test of 
the predictions in Ref. 45. Transmission and reflection 
MIR magneto-optic meaurements on Gai-^^Mn^^As ran- 
dom alloy fllms in Figs. 10 and 11, respectively have 
been recently measured in this MIR range and support 
the presence of such strong magneto-optic signal in the 
MIR regime. -^^^ 
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FIG. 9: Faraday and Kerr effects for x — 6% Mn concentra- 
tion and p = 0.4 nm"^. After Ref. 57 



Figure 10 plots the real a) and imagnary b) parts of 
the MIR [9p] probed at a photon energy of 118 meV and 
a sample temperature of 7 K. The hysteresis is weak and 
similar to out-of-plane dc magnetization measurements 
performed on this sample (a 6 % Gai_a;Mn2^As random 
alloy film). Both the real and imaginary parts of [9y] 
appear to be proportional to the film magnetization, as 
expected. The Rc[0f] signal at magnetization saturation 
calculated by Sinova et al.'^^ are on the order of 0.04 rad 
for a film with a similar thickness and Mn concentration. 
However such calculation are done in optimally annealed 
samples and further experimental measurements in this 
more metallic regime will be helpful to establish the va- 
lidity of the model. 

In Fig. 11, a two-step hysteresis loop is observed both 
the MIR (118 meV) Re[6'K] as well as the dc magneti- 
zation for a 1 % Gai_j;Mna;As random alloy film. This 
behavior is similar to that observed in MCD at a photon 
energy of 2.83 eV by Ref. 108. The authors in Ref. 108 
claim that such behavior cannot be attributed to GaMn 
clusters, and therefore is a signature of the Gai_a;Mnj;As 
alloy. 

The MIR magneto-optical measurement technique 
used in the above measurements is described in detail 
in Ref. 123, and results using this technique on non- 
magnetic systems are presented in Refs. 124 and 125. 
The results shown in Fig. 10 and 11 serve to demon- 



FIG. 10: RelOp] a) and Imlep] in the MIR range b) as a 
function of applied magnetic field for GaAs(6 % Mn) random 
alloy epilayer at 7 K and probe radiation energy 118 meV, 
after Ref. 122 



strate the experminetal feasibility of such studies, but 
further work is required to interpret these results and to 
compare them with theoretical predictions. 



V. SUMMARY 

This review has focused on highliting the rapid devel- 
opments in material research of metallic ferromagnetic 
(III,Mn)V semiconductors over the past few years. Al- 
though the successes have been extensive in the under- 
standing and development of the bulk properties of the 
(Ga,Mn)As and (In,Mn)As systems within the metallic 
regime, many challenges remain both in other (III,Mn)V 
based semiconductors and lower dimensionalities such as 
digitally doped heterostructures. Both transport and 
magneto-optical effects will remain a very fruitfuil ground 
to make progress in these materials. At the time of writ- 
ing, there is still a lot of work being done trying to sort 
out the more complicated physics of (Ga,Mn)N and sev- 
eral reports of high Tc measurements in the digitally 
doped materials. We simply keep watching in fascina- 
tion the rapid developments and challenges posed by new 
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tion (courtesy of S. Wang) for GaAs(l % Mn) random alloy 
epilayer. After Ref. 122 



experiments which will undoubtedly generate interesting 
physics and new understandings of these fascinating ma- 
terials. 
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